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Spontaneous four-wave mixing can generate highly correlated photon pairs from atomic vapors. 
We show that muhi-photon pumping of dipole-forbidden transitions in a recoil-free geometry can 
result in ultra-bright pair-emission in the full An solid angle, while strongly suppresses background 
Rayleigh scattering and associated atomic heating. Such a system can produce photon pairs at rates 
of ~ 10^^ per second, given only moderate optical depths of 10 ~ 100, or alternatively, the system 
can generate paired photons with sub-natural bandwidths at lower production rates. We derive a 
rate-equation based theory of the collective atomic population and coherence dynamics, and present 
numerical simulations for a toy model, as well as realistic model systems based on '^^^Cs and ^^^Yb 
level structures. Lastly, we demonstrate that dark-state adiabatic following (EIT) and/or timescale 
hierarchy protects the paired photons from reabsorption as they propagate through an optically 
thick sample. 

PACS numbers: 42.65.Lm, 42.50.Ar, 42.50.Dv 



I. INTRODUCTION 

The study of correlated/entangled photon paks has 
long been a central topic in the field of quantum optics 
[T]. The importance of paired photons is two-fold: they 
i) provide powerful tools to test the peculiar aspects of 
quantum mechanics, such as violations of local-realism 
[21 [31 13]; and ii) they hold promises for advancements 
in quantum measurement, communication, and informa- 
tion processing H [THH [3 [THl [II]. Over the past 
few decades, spontaneous parametric down-conversion 
(SPDC) in nonlinear crystals has been the standard 
source of photon pairs [HI |T3|. More recently, an al- 
ternative class of biphoton sources has emerged, based 
on optical four-wave mixing (FWM) in atomic vapors 
[Tl[Til[I51[I71[TSl[Tni[2ni[2Illia- These approaches 
rely on collective effects [231 [23] to greatly increase the 
probability of correlated emission events. Compared to 
SPDC, photon pairs generated via FWM in general have 
a much narrower bandwidth, significantly greater tempo- 
ral and spatial coherence, and much higher conversion ef- 
ficiencies. They are thus particularly suitable for hybrid 
quantum communications and computations employing 
atoms and photons |lll 125] , and for high-precision quan- 
tum measurements and imaging [HJ [TU] . 

At present, FWM photon pair sources can be catego- 
rized into three types by level configuration. The first 
type, built on atomic two-level systems, is a connected 
double-Rayleigh emission process [26] [27] . Due to strong 
background Rayleigh scattering, however, the resulting 
pair correlation is very weak, without satisfying the nec- 
essary Cauchy criteria for biphoton correlation [IS] [H] . 
A second type is configured on two-photon cascade emis- 
sion in a four- level system [50] [2H]- While high-fidelity 
photon pairs are generated, due to the unequal wave- 
Icngthcs of two cascade photons, the phase- matching con- 
dition for collective emission can only be satisfied if the 
first photon is emitted by chance into a specific small 



solid-angle, thus unpaired emission dominates the over- 
all radiation, resulting in a relatively low conversion effi- 
ciency. The third type employs Raman FWM (hereafter 
referred to as "RFWM") in multilevel systems, config- 
ured on double-A [H [TTl [TSl [H] [22] or "X" ^ level 
diagrams. The major challenge in these schemes is to 
suppress background Rayleigh scattering, which tends to 
rapidly overwhelm paired emission. Three approaches 
have been proposed to for this suppression, including i) 
using frequency selectors to filter out Rayleigh photons 
[TT] : ii) collecting pairs along emission directions where 
the dipole pattern leads to zero Rayleigh emission |18j : 
and iii) using a single-mode optical cavity to suppress 
Rayleigh transitions [Hj. While yielding up to 10^ pairs 
per second, all of these setups are unidirectional, where 
photon pairs are produced only along certain directions. 
This restricts the obtainable beam brightness of the pho- 
ton pairs, since in each momentum mode, the time sep- 
aration between pairs must be large relative to the cor- 
relation time. Lastly, in aforementioned FWM schemes 
where unpaired emissions dominate, atomic samples are 
rapidly thermalized due to random atomic recoils, limit- 
ing applications of these schemes to 'hot' vapors only. 

Background Rayleigh scattering occurs when there is 
a spontaneous one-photon channel by which atoms can 
return to the initial internal hyperfine level without com- 
pleting the desired biphoton emission cycle. We propose 
to eliminate this in RFWM by replacing the single pump 
laser with a multi-photon pump process that drives a 
one-photon dipole-forbidden transition. The second part 
of our proposal is to use 'recoil-free' pumping, meaning 
that the k-vectors of all of the driving fields sum to zero. 
Phase matching, which is enhanced by collective effects, 
will then be satisfied whenever the two paired photons 
have equal and opposite momenta, allowing paired emis- 
sion into the full 47r solid angle. As a result, the biphoton 
emission rate will be enhanced by up to four orders of 
magnitude over unidirectional approaches. Furthermore, 
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since Rayleigh scattering has been eliminated, atomic 
thermalization will be strongly suppressed. Thus this 
scheme can also be applied in ultracold vapors including 
Bose-Einstein condensates (BECs), for which it might 
also be viewed as a novel in-situ non-demolition imaging 
technique. 

A brief organization of this paper is as follows. In sec- 
tion [n] we study a 'butterfly' biphoton protocol using a 
simplified model, which is a viable simplification of real- 
istic configurations. In section |III[ we give two realistic 
implementations, one employing the 399nm-line transi- 



tion in Ytterbium (section III A ) , and the other using 



the 852nm D2-line transition in Cesium (section 
In section IIVI we then address the critical issue of 



IIIBl 



how to 

avoid reabsorption as the photon pairs propagate in an 
optically thick sample. This is followed by a discussion 
and conclusions in section IVl 



II. BASIC MODEL 

In this section, we will use a simplified 'toy model' 
to show the important physics of the butterfiy scheme 
for ultra-bright photon-pairs. We present the schematic 
model in section |II A[ and employ a set of rate equations 
to solve for atomic dynamics in section [lIB[ Then, we ex- 
amine the time and polarization correlation of generated 
photon pairs in |II C| and |IID[ In section |IIE| we esti- 
mate the threshold temperature of the present biphoton 



A. Toy Level Scheme 

A schematic level diagram of the butterfiy scheme is 
shown in Fig. [l](a). While greatly simplified with re- 
spect to a realistic level-scheme, this model will serve to 
illustrate the important dynamical effects. The physical 
mechanism is best illustrated in a quantum trajectory 
picture, as follows. A sample of atoms initially pumped 
into the |1) state, is first weakly coupled to the excited 
|2) level via a multi-photon pump process. This imparts 
a net recoil momentum of /iK, so that for an initial 
momentum fi.q, an atom excited to |2) has a momen- 
tum of /i(q -I- K). This excited atom will then sponta- 
neously decay to |3), emitting a 'signal' photon with a 
random momentum Kk, shifting the atom's momentum 
to h{q + K — k). Decay from |2) back to state which 
would generate background Rayleigh scattering, is for- 
bidden by dipole selection-rules. The atom in state |3) is 
then rapidly repumped to |4) by a strong multi- photon 
coupling process. The coupling fields are arranged to 
yield a net momentum of —UK., leading to a momentum 
of /i(q — k) for the atom. From |4), the atom decays back 
to the |1) state, emitting an 'idler' photon. Collective 
effects will strongly enhance the scattering probability if 
the atom can be returned to its initial momentum state 



hq, which will result in the idler photon being emitted 
with momentum — fi,k. 

The collective enhancement mechanism can be under- 
stood by noting first that the emission of the signal pho- 
ton with momentum hk. imprints 'which atom' informa- 
tion onto the atomic ensemble via atomic recoil, provided 
of course that the single-photon recoil momentum is 
larger than the momentum coherence length of the sam- 
ple (which for a thermal gas of free particles is the inverse 
sample length). If the idler photon is then emitted with 
phase-matched momentum, hki ~ — fik, then the atom is 
restored to its initial momentum state of fi.q, thus 'eras- 
ing' the 'which atom' information, so that many-body 
interference enhances the emission rate by a factor N. 
The 'which atom' information will be effectively erased 
only when the magnitude of the momentum difference be- 
tween the initial and final states is less than or equal to 
the momentum coherence length of the atom, |k + k^j < 
kcoh, which leads to a collective emission solid-angle of 
^coh ~ '^{kcoh/k)'^ ■ For a lone atom, and neglecting 
the dipole emission pattern, the probability of correlated 
emission, i.e. the probability of kj falling within flcoh of 
— k, is /'c(l) = f^coh/47r with a non-correlated probabil- 
ity P„c(l) = 1 — ^coh/4:TT. For a sample of N atoms, the 
differential probability inside ^Icoh is enhanced by A'' rel- 
ative to the differential probability outside ^coh, leading 
to a collectively enhanced correlated emission probability 
ofP,{N) = Nn,,h/{Nn,oh+'^Tr-n,,h). The probability 
of non-correlated emission is correspondingly reduced by 
unitarity to P„c(A) = {An - Vtcoh) / [N^coh + 47r - Vtcoh)- 
Under the condition N^c.oh/^T^ ^ 1, the probabilities are 
then given approximately by Pc{N) = 1 — Air/N^lcoh ^ 1 
and Pnc{N) « Air/Nilcoh ^ 1- For a thermal sample 
of dimension L, we have kcoh ~ ^/L, so that flcoh 
7r/(fci)^. This gives a correlated emission probability of 
Pc = l-1/D, where D = N/{2kLf = nA^L/lGTr^ is the 
optical depth of the sample. Other considerations aside, 
we clearly see that strong photon pair-correlations in k- 
space are achieved by maximizing the optical thickness 
of the sample. 

The fact that the driving and coupling fields have zero 
net momenta makes the scheme in a limited sense 'recoil- 
free', so that phase-matched collective emission can occur 
regardless of which direction the signal photon randomly 
'chooses'. This is illustrated in Fig. [l](b), where for an 
arbitrarily chosen k, the atomic dynamics in the space of 
recoil momentum undergoes a closed, diamond-like cy- 
cle. This distinguishes the present setup from competing 
biphoton protocols employing non-counter-propagating 
driving and coupling, where pair emissions are restricted 
to be in a plane perpendicular to the nonzero net mo- 
mentum [IillTHllMlliH]. 



B. Population Dynamics 

To study the system's dynamics, we quantize the 
atomic center-of-mass motion, onto the eigenmodes. 
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FIG. 1; (Color online) A schematic model of the butterfly 
scheme. Figure (a) draws the simplifled level diagram, em- 
ploying multi-photon driving and coupling pumps. The no- 
tation I j, k) indicates a single atom state in internal level | j) 
and with momentum Kk. The momentum fiK is the net mo- 
mentum of the multiphoton transition from levels |1) to |2). 
Figure (b) shows how a phase-matched diamond-like cycle 
in atom-recoil momentum space exists for any signal photon 
emission direction. 



{|q)}, of a box of dimension L, where L is the ensemble 
dimension. As the box traverse time, L/v, is long com- 
pared to the relevant dynamic timescale, (r)~^, we can 
safely impose periodic boundary conditions, so that the 
allowed q values lie on a three-dimensional cubic lattice 
with spacing 27r/L. Eliminating the scattered light-field 
via the Markoff approximation, and taking into account 
the exchange-type symmetry respected by the Hamilto- 
nian and the initial conditions, allows us to derive a set of 
rate equations for the atomic population dynamics. The 
advantage of rate equations, as opposed to a mean-field 
approach, is that spontaneous decay is incorporated via 
the usual (A^ + l) factors, thus avoiding the need for noise 
operators. For the toy model depicted in Fig. 1 (a), the 
rate equations are: 

^iVi = {^^921 - c.c) + ^r4/3q4A^q4(A^l + 1), (1) 

q 

= (^!d^^21 - c.c) - ^r2/3q2iV2(iVq3 + 1), (2) 

q 

+ ^621^ [ri/J,4Af,4 - rjf),j(«,j + 1)1 , (3) 

q 

^iVq3 - {n,g^43 - c.c) + r2/3q2iV2(iVq3 + 1), (4) 
d i 

j^N^4 = ^2 (^^^q43 - c.c) - TiPqiN^iiNi + 1), 

- r4(l - /3q4)iVq4 (5) 

d .^c,^, ^ 1 

[r2/3q2A^2 - r4/3q4(iVi + 1) - r4(i - p^i)] ,(6) 



a detailed derivation of these equations, as well as pre- 
cise definitions of the variables, is presented in appendix 
\K\ Oversimplifying slighty, we can think of Ni and N2 
as the populations of state |1) and |2), respectively, with 
g2i being the corresponding coherence operator. iVq3 is 
the expectation number of atoms collectively excited into 
state 1 3) via emission of a signal photon with k « — q. 
Similarly, 7Vq4 is the expectation number of these atoms 
transferred to state |4) by the coupling laser, while 9^43 is 
the coherence between these two collective states. Lastly, 
and fie are the effective Rabi frequencies of the driv- 
ing and coupling transitions. The spontaneous emission 
rates for |2) — > |3) and |4) |1) decays, are r2 and 
r4, respectively, and /3q2 and /?q4 are the branching ra- 
tios for emission into the coherent-emission solid angle 
^coh = T^/{kL)'^, with respect to the ±q directions. For 
a spherical sample of radius L/2, a careful calculation 
p9] gives 

/3qM = (l-|'7-4n|^(^)'; M = 2,4; (7) 

where d2, and d^ are the unit vectors along the dipole 
moments of the |2) — > |3) and |4) |1) transitions. Aside 
from dipole-emission factor, (1 — \q- dj\'^), this is in good 
agreement of our initial estimate ric/47r ~ /{AttL)'^. 
The quantity £'q4 — NiPq^ is then the optical depth 
along the ±q directions, with respect to the |1) |4) 
transition. 

Strictly speaking, in the toy model, atoms which spon- 
taneously decay from |4) to |1) by emitting rogue photons 
can still participate in the next-round pair-emission cycle, 
as the initial momentum of the atoms in state |1) is irrele- 
vant to the collectivity (in the Doppler-free regime). This 
recycling process can be included in the rate-equation 
model by inserting a re-feeding term X)q 1^4(1 ~ /3q4)-^q4 
to equation (ll]). In this case, the total atom number 
iVi -I- iV2 + JZq -^q3 + -^q4 ^i^^ conserved, while still 
correctly describing the emission of rogue (uncorrelated) 
photons. In realistic schemes, however, atoms can also 
spontaneously decay to other ground levels which are not 
shown in figure [l] and/or they can decay from interme- 
diate pumping levels to |3). These atoms will not be 
able to participate in further collective emission cycles, 
unless the are somehow repumped back to |1). Thus to 
avoid overly optimistic predictions, we have chosen not 
include the recycling process in the rate-equation model. 
It is noted that by completely excluding the recycling 
process, the pair generation rate is underestimated, as a 
fraction of rogue-photon emitting atoms will always be 
recycled. On the other hand, atoms which decay from |4) 
into levels other than |1) may be repeatedly re-excited to 
other levels by the pumping fields, and thus emit addi- 
tional rogue photons. If such rogue photons can not be 
filtered out, they will contribute to the impurity of the 
collected biphoton beams. For present, however, we only 
focus on the short-time behavior of the system, up to a 
point when the atom loss is about 10%. In this time in- 
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terval, the present dynamical model with rate equations 
([T])-([6]) is reasonably valid. 

The total emission rates for signal and idler photons, 
corresponding to the (enhanced) decay rates of the |2) 
and 1 4) levels, are 

Rs - r2^/3q2iV2(iVq3+l), (8) 

q 

Rl = r4^/3q4iVq4(^l + l), (9) 

q 

respectively. Assuming steady-state, clearly we must 
have Ri < Rs- If Ri < Rs, more signal photons are 
generated than idler photons, so that pairing is weak. 
Thus, at a minimum, strong pairing requires Ri — Rs- 
Focusing on a single q mode, and assuming /3q2 — /3q4 
and r2 = we find 



Rijk) 
RsCk) 



Ni + 1 
No 



N, 



q4 



N. 



q3 



1 



(10) 



If we assume a strong drive, fi^ ^ r2, we have N2 ^ Ni, 
which means we must also have A'q4 ~ iVqs ^ 1, which 
in turn requires a strong coupling field, flc ^ TiPqiNi- 
We find, however, that dynamically this approach doesn't 
work, as it leads to a build-up of population in A^q3 and 
A'q4 without strong pairing. This leaves the case of weak 
driving, ild ^ r2 so that N2 <^ Ni- This then requires 
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Ni 



the condensate only depletes at a rate of Rrogue, so that, 
e.g. A^_D/10 ~ lOiV photon pairs could be generated 
with only 10% of the condensate atoms being lost. This 
means that if desired, the present butterfly scheme can be 
used to directly image condensates in situ, in a relatively 
nondestructive manner. In other words, the BEC would 
exhibit resonance fluorescence, but with strongly sup- 
pressed heating. The additional brightness might, e.g., 
yield improved atom-number estimation. 

Similarly to the case of two-photon cascade emission 
[5] , here the bandwidth of both signal and idler photons 
in the strong coupling regime (|f2c| 3> DTn) is given 
by r4D/2, where Y4^D is the superradiance broadened 
linewidth of state |4). This is because the intermediate 
1 3) and |4) levels participate the pair emission by first 
forming dressed states of |±) — (|3) ± |4))/'\/2. The pro- 
cess of pair emission depicted in figure[l]is then physically 
equivalent to cascade emission from |2) to |1), through ei- 
ther |-|-) or |— ) states. In further analogy to cascade emis- 
sion, the frequency-sum of the signal and idler photons 
has a much narrower linewidth of |f2dp/r2, as given by 
the reciprocal of the pair emission time. In the weak cou- 
pling regime, \Q,c\ <C DT4^, the coupling dynamics is over- 
damped, so that the |4) level can be adiabatically elimi- 
nated. Then the pair emission is effectively through the 
1 3) level only, the linewidth of which is broadened from 
zero to |ricP/-Dr4 by the coupling laser. This, together 
with the fact that we choose \Q.c\'^ / DT^ » |rJdP/r2 for a 



which can be arranged by adjusting strong pairing effect, gives a bandwidth of ^fJ- < D r4 



the drive and coupler strengths and detunings, and pro- 
vided A'^q4 <C 1, results in strong pairing [TS1[22]. 

The impurity of collected biphoton beams in the 
present toy-model comes from spontaneous emission of 
rogue idler photons into non phase-matched angles. 
The total emission rate of rogue photon is given by 

Rrogue = r4l]q^q4(l " Pt^i) « ^iY.c^^^■i^ where 

/3q4 <C 1. The ratio of paired idler to rogue idler pho- 
tons is then given by Ri/ Rrogue = {Ni + l)/3, where 
P = X]q/3q4^q4/X]q^q4 is the mean collectivity aver- 
aging over emission angles. For a spherical cloud of ra- 
dius R, this is roughly X^/IGttR^- Defining the optical 
depth of a spherical cloud a,s D = Ni(3, the pair to rogue 
ratio is simply Rj/ Rrogue — D, i.e., the optical depth. 
For typical samples of D '--^ 100, there is then about one 
rogue photon per 100 pairs. 

In general, the emission of a rogue photon leads to 
heating of the sample, due to random non-zero net-recoil. 
In a unidirectional scheme, the vast majority of photons 
are not paired, so that heating occurs at the usual single- 
photon decay spontaneous heating rate. Such schemes 
are thus only applicable to samples well above the re- 
coil temperature. With the omnidirectional approach, 
the majority of spontaneous photons come in correlated 
pairs and thus impart no recoil kick. The heating rate is 
then reduced by a factor of the optical depth D, which 
should allow interesting experiments to be performed at 
or below the recoil temperature. For example, if a BEC of 
iV atoms is used as an omnidirectional biphoton source. 



for the emitted photons. In the very weak coupling limit, 
with \^c\ < Vd^a, this means that correlated photon 
pairs of sub-natural bandwidth can be generated. 

To study the performance of this system, we numeri- 
cally solve equations ([l])-([6]). One question that must be 
decided is how-many and which box eigenmodes need to 
be included in the simulation. Firstly, we note that in 
momentum space, the width of the energy shell for the 
signal photons is Ak = T4D/C ~ lm~^, while the size of a 
single eigenmode is l/L- Thus for normal sample sizes of 
L <C Im, only a single shell of modes of radius 27r/A par- 
ticipates in the dynamics, where A ~ 10~''m is the wave- 
length of the |2) |3) transition. This leads to a mode 
number of M = Att/Uc ~ (L/A)^, which could range any- 
where from 10^ to lO^'', depending on the sample size. To 
handle the dipole pattern, we then sort these modes {q} 
into 15 groups. The j-th {j — 1,2, ...15) group contains 
those quasi-modes satisfying 9q E [Sj,Oj + 7r/30), where 
9j = {j — l)7r/30, and ^q e (0,7r/2] is the angle between 
q and the d2 axis. Modes in the j*'' group, are assigned 
with mean collectivity parameters, fj2, and fj^, obtained 
by averaging over modes inside the group. In this ap- 
proximation, all quasi modes within the same group will 
yield identical dynamics, so that only of one mode in each 
group needs to be included in the dynamical model. The 
effects of the other modes can then be included by weight- 
ing each representative mode by the number of modes on 
the interval [0j - 7r/30,6'j) ~ (i/A)^ sin(j7r/30). In this 
way, we are able to reduce the number of coupled rate 
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equations from 10^-10^", down to fewer than 100, while 
still incorporating the effects of the dipole radiation pat- 
tern. 

We consider an example of = 10^ atom, with a 
spherically-symmetric Gaussian density distribution of 
radius L = 34A, corresponding to an optical depth of 
D = 70. We consider resonant driving and coupling 
pumps propagating along z and — z directions, and take 
T2=Ti = r, fid = O.ir, and = lOOF. The results 
are shown in Fig. [2] (a)-(d). In figure (a), we plot the 
time evolution of Ni,N2, showing that at all times only a 
small fraction of atoms are excited, i.e. N2 ^ Ni, as re- 
quired. In figure (b), we plot the mean atomic collective 
excitation numbers and obtained by averaging 
over q. Both are found to be of order of 0.01, so that 
there is negligible overlap between subsequent pairs in a 
given mode. The total photon-pair number and lost atom 
number are shown in figure (c), where they are found to 
increase linearly in time with fitted rates of 8.3 x lO'^F 
and 1.4 X lO^r. 
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FIG. 2: (Color online) Figures (a)-(b) show the evolutions of 
Ni, N2, N3, N4,, while (c) compares the number of generated 
photon pairs Npair and lost atoms Nioss ■ Figure (d) plots the 
second-order correlation function g''^-' (k, — k, r) for k = fci 
and ]i. — kx, respectively. Parameters are given in text. 



C. Photon Pair Correlation 

To estimate the time correlation of the two photons, 
we calculate the time-averaged second-order correlation 
function g^^)(k, — k, t), defined as j2] 

g^'\k~k,r)^^ £ dtx 



(aL(*)ak.(t))(a^_k»(i + r)a^^^^{t + t)) 



,(11) 



where Oks and are the annihilation operators for sig- 
nal and idler photons, which may or may not differ in po- 
larization for a given k, and T is the averaging window. 
The correlation function is evaluated by first using adia- 
batic following to write the photon operators in terms of 
the atomic operators, where 



ak.(t) = ak.(0)e-*"''* - igk 



sm[{LOs - ^k)V2] 



xe-*("=+">')*/2a2.k3, 



aUt) = ak.(0)e-*">'* - Zffk 



sin[(ws - Wk)V2] 
{ujs - Wk)t/2 



xe-*('^=+'^'')*/2CTk4,i- 



(12) 



(13) 



Here, gk is the atom-photon coupling constant, Wk ~ 
c|k| is the photon's frequency, and is the resonance 
transition frequency. The collective atomic operators are 
defined as 

I72.k3 = ^3,2 ^^f{^ + ^- Q)'S'Q2q3 (14) 

q,Q 

0-k4,l = 0'I,k4 = Q)'^Q4ql (15) 

q,Q 



where the ensemble operator S^i, and the structure func- 
tion /(k -t- q — Q) for the atomic sample are precisely 
defined in appendix [A| Eq. (AI8I and (A20l, respec- 
tively. Inserting this result to the correlation function 



(111 gives 



ff(2)(k,-k,r) 



1 
T 



dt X 



(0'k3,2(^)'7l,k4(^ + T)o-k4,l(^ + 'r)q-2,k3(^)) 
('5'k3,2(i)<5'2,k3(i))('5'l,k4(i + T)(Tk4,l(< + t)) 



,(16) 



where any contribution from vacuum electromagnetic 
(EM) fiuctuations, given by aks(O) and aks(O), vanishes 
when tracing over the EM vacuum. 

As demonstrated in detail in Appendix|X] because our 
system is in a symmetric collective state throughout the 
dynamics, the products of coherence operators can be 
written as the products of number operators. Applying 
this to equation ( 16 1 , we have 



J2) 



dt 



(gq3,2(^)iVq4(^ + T)a2,q3(t)) 

N^i{t + T)N2{t){N^3{t) + \y 



(17) 



where we have replaced the number operator Ni with its 
meanfield A^i and approximated N^ + l « A^i . This result 
is a valid approximation as the |1) level is macroscopically 
occupied, with iVi ^ 1. The discrete momentum q is the 
nearest neighbor of k, with which |k — q| is minimized. 
Following the standard procedure [2j , we write 

2 4 

+ [vi{r)gi2{t) + m{r)gci34{t) + c.c] , (is) 
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where we have used the fact that the coherent terms 
(Pq32) and (^q4i) are zero throughout the dynamics, since 
neither pumping nor purely spontaneous/superradiant 
decay processes will generate such coherences. The coef- 
ficients Xj(T) and riji^r) are determined by studying the 
linear response of Nq4{t + r) to small perturbations in 
each variable, via numerically solving the rate equations 
([T])-([6]). For instance, to determine X3(t), we apply a 
small, instant perturbation i5A^q3 to the system dynam- 
ics, i.e., making iVq3(i) Nq3{t) + SNq3. We then calcu- 
late the resulting drift SNq^it + r) from the unperturbed 
value 7Vq4(i + T), with which the coefficient is determined 

as X3(t) = '- 

FoUowin 



the quantum 



regression theorem [30^, in 
equation ( 17 1 we write Nq^^t + r) at time t + t in terms 
of operators at time t, using the result (18 1. The re- 



sulting expression is then normalized and factorized into 
products of the occupation numbers, giving 
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(2) 



1 



f Xajr) ~ X2(t ) , rt{(T)Q4i(t) 
V ^q4(i + T) 



'72(T)e'q43(t) 



(iVq3(i) + l)A^q4(i- 



iV2WAfq4(i + r) 
(19) 



As seen in figure [2] (a) and (b) , the atomic dynamics un- 
dergoes quasi-steady-state evolution, with slow damping 
due to atom losses. This allows us to approximate the 
populations and coherences as constants for the time pe- 
riod of interests, while effectively taking T ^ oo. 

As seen in equation (19), a strong time correlation be- 



tween signal and idler photons requires 7Vq4 <C 1. Phys- 
ically, this is because the one-to-one correspondence be- 
tween signal and idler photons will be spoiled if there is 
more than one 'atom' in the same collective mode at the 
same time. This result agrees with the established ex- 
perimental criteria, where the driving field must be weak 
such that most atoms remain in the initial ground state 
|17l [18J. In this weak driving regime, we find to a good 
approximation that 



a 



(2) 



1 



1 



< ^q4 > 



(20) 



with < A'q4 > being the time-averaging value of N^^. 
For a weak coupling between |3) and |4) with fie ^ 



« 1, 



which leads a 
A 



£)r4, we find xzij) ~ u^.^ 

time delay iff-p- between signal and idler photons 

strong signal-idler correlation with g'^^ ^ 1 then re- 
quires iVq4 to be very small. Since the pair generation 
rate i?/ '-^ -^q4, this would ultimately limit the achiev- 
able beam brightness of photon pairs. For a strong 
coupling with fie DV^, however, we find X3(t) ~ 
sin^ (j^ct) exp(-i£)r4T) 1. In this relatively 
larger A^q4 can give the same correlation, thus elevating 
the obtainable beam brightness. The correlation function 
in this case exhibits oscillatory and damped behaviors 
with sharp peaks, associated with Rabi-oscillations be- 
tween 1 3) ^ 1 4) [17]. The time delay between signal and 



idler photons is roughly (D r4)~-'^. We note in all regimes, 
the bandwidth of signal and idler photons is the recipro- 
cal of the delay time. In figure |2](d), we plot the second- 
order correlation functions for photon pairs propagating 
along ±z and ii directions. Both cases exhibit sharp 
peaks of widths ~ O.OSF^^, indicating strong temporal 
correlation which violates the standard Cauchy-Schwartz 
inequality by a factor > 1000. 



D. Polarization entanglement 

We now examine the polarization entanglement of 
paired photons. We consider a pair of signal and idler 
photon individually travelling along k and — fc directions, 
with k{9,4>) = sin 9 cos (px + sin 9 cos (py + cos 9 z , where 
z is along the K direction. For ^2,4 = zk iy), the 

probabilities for signal photons to be left and right cir- 
cularly polarized along k, denoted as and en, are 
/3^{9) ^ (l-hcot^l)"^ and f3§{9) ^ /3^{tt ^ 9), re- 
spectively. Similarly, for the idler photons we have 
/3f = /3|(6»), and /Jf = /3f (6I). The probabihty for pho- 
tons to be in opposite circular polarizations along k (thus 
in the same polarizations along each's propagating direc- 
tion) is then 



P{9) = P^{9)l3^{9)+l3^{9)/3f{9) = 



1 



cot'' I 



(l-t-cot^f) 



(21) 



As seen in Fig. [3] this is extremely fiat around 9 — 0, tt, 
where P{9) w 1 — |mod(0, tt)^ w 1, meaning that photon 
pairs emitted over a wide range of 9 will yield strong 
polarization entanglement. Due to the temporal overlap 
of signal and idler photons, each pair emitted within the 
strong correlation angle is approximately in the Bell state 
of = ^(le^ei) -I- Iciefl)). For example, pairs with 
one photon emitted within 9 < 0.5 (corresponding to 
17% of total emitted pairs) have an entanglement fidelity 
> 99%. 




FIG. 3: The probability P{6) for paired pliotons to be in 
opposite polarizations as a function of emission angle 9, cal- 
culated for the toy model. 
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E. Critical Temperature 

A biphoton source in an atomic ensemble relies on 
collective enhancement (superradiance) to increase the 
probability of correlation emission events, i.e. enforce 
'phase matching' for the light waves. A necessary 
condition for collectivity is that the Doppler linewidth 
is smaller than the photon's superradiance-broadened 
linewidth, the latter of which is a factor D greater than 
the natural linewidth [3I| 132] . Physically, this condition 
is required to suppress dephasing of the collective atomic 
excitation during the superradiant emission process. The 
present scheme is recoil-free in the sense that the driv- 
ing and coupling fields have net zero momenta, so that 
phase matching requires only that the signal and idler 
photons have equal and opposite momenta. However, an 
atom undergoing the biphoton emission cycle does pos- 
sess a nonzero recoil momentum while in the intermediate 
states |2), |3) and |4). Thus the system is indeed subject 
to some Doppler broadening , as opposed to conventional 
Doppler-free geometries where true two-photon pumping 
creates an intermediate excited level with zero recoil mo- 
mentum |33j . 

A useful alternative picture of Doppler broadening has 
emerged from considering superradiant process in ultra- 
cold atomic gases [34l|35], based on the interplay between 
the spatial coherence length and the recoil velocity. In 
this picture, a thermal atom is viewed as a spatial 'blob' 
of coherence, whose size is given by the thermal coher- 
ence length Xcoh = h/^j2m,kBT. When a single photon is 
collectively absorbed, each atom is placed in a quantum 
superposition of its initial state, and the excited state, 
which is also a coherent blob, but one which is moving 
at the recoil velocity relative to the initial blob. If col- 
lective excitation lives too long, then the initial and ex- 
cited coherence blobs no longer overlap in space, at which 
point it no longer matters whether or not they overlap 
in momentum space. Overlap in both momentum and 
position space is required in order for the 'which-atom' 
information to be erased. This sets a minimum criterion 
for collective effects as, VrTc ^ Xcoh, where — hK/M 
is the recoil velocity, and Tc is lifetime of the collective 
excitation. With = I/t^, this gives 3> 7d{T) as the 
necessary condition for collective enhancement, where 



Td{T) = K 



(22) 



is the usual Doppler broadening contribution the 
linewidth. We therefore see that Doppler broadening 
becomes significant only when Td{T) is greater than or 
comparable with the linewidth of signal and idler pho- 
tons, both given by D The threshold temperature of 
the system is then determined as 



(23) 



10^s-\ 10h-\ and lO^s-i, we have = 10, 10^, and 
10* kelvin, respectively. Due to the scaling of Tc as Ff, we 
see that there could be a significant advantage to using 
an atomic transition with a large natural linewidth. 

This means the present butterfly scheme can be imple- 
mented with room-temperature atomic clouds, excluding 
Cs which has F ^ 10^, and so must be cooled below 
lOK. We note that at i ntermedia te temperature, defined 
as DF4 > Td{T) > + r|, where A is the detun- 

ing of ild (which so far is taken to be zero), the effective 
linewidth of level |2) is broadened from F2 to Fd(T). For 
A = 0, this is the case when T > 1 kelvin. The only 
consequences of this, however, are that the various laser 
intensities and detunings would have to be adjusted ac- 
cordingly to maintain the condition (10 1, and the sum- 
frequency linewidth would be similarly Doppler broad- 
ened. The pair correlations, however, will be unaffected. 
In contrast, for T > Tc the pairing effect will completely 
disappear. 



III. REALISTIC MODELS 

In the above section, we used a simplified toy model 
to illustrate the physics of the butterfly scheme for ultra- 
bright photon pairs. In this section, we provide two re- 
alizations of the butterfly scheme, implemented with Yt- 
terbium (section III A I and Cesium atoms (section III B I , 
respectively. The main difference between toy and real- 
istic models is that in the toy model, we use effective 
Rabi frequencies to describe the multi-photon pumping 
dynamics. In practice, this is an over-simpliflcation of a 
complex dynamical process, and excludes the background 
Rayleigh/Raman scatterings associated with the inter- 
mediate states. To study these features, in the following 
section, we shall expand the toy dynamical model to ac- 
count fully for the multi-photon nature of the drive and 
control fields. 



above which collective effects disappear. Taking typical 
|K| = lO^m-i, D = 100 and m = lO-^^kg, for F4 = 



A. Ytterbium Atoms 

We now consider a realistic butterfly level scheme con- 
figured on the 399nm-line of the 6s^ ^So<-^ 6s6p ^Pi tran- 
sition in ^'^-'^Yb atoms, as shown in figure H (a). The 
atoms are prepared in the |1) = \F = 1/2, mp — —1/2) 
state and then follow a FWM cycle which deposits them 
in level |3). A second independent FWM cycle then re- 
turns them to the initial |1) state. The driving FWM 
cycle consists of one violet laser and two infrared lasers, 
with Rabi frequencies 51 1, and fis. The coupling 
FWM cycle similarly contains one violet and two in- 
frared lasers, with Rabi frequencies of fl^, fis and fig- 
The net momenta of the driving lasers are equal and 
opposite with that of coupling lasers. Signal photons 
are emitted as |2) = |F = 3/2, nip = 3/2) atoms 
spontaneously decay to the only dipole-allowed state of 
1 3) = \F = 1/2, mi? = 1/2), and idler photons are gen- 
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erated as each atom in |4) = \F = 3/2, mp = —3/2) de- 
cays collectively back to |1) with initial momentum /iko. 
Both transitions have a natural linewidth of 0.2 GHz. 
Atoms which spontaneously decay to other momentum 
modes of level |1) generate unpaired rogue photons. Yet, 
they are still able to participate in the next-round col- 
lective emission cycle. In contrast, atoms which sponta- 
neously decay from the intermediate pump levels to |3) 
are not eligible for the collective emission. They nonethe- 
less can be repumped back to |1) in subsequent dynam- 
ics, during which more rogue photons will be emitted. 
For simplicity, however, in the present model we treat 
all unpaired-emission events as permanent atom losses, 
thereby excluding any of the aforementioned recycling 
or repumping processes. We note here the background 
scatterings, involving spontaneous single-photon decay of 
mp = — hyperfine states of the 6s6p level, are 
efficiently suppressed by making fli and fl^ far-detuned 
from resonance. The intermediate 6s5d level has a long 
lifetime of 6700ns, compared to 5ns for the 6s6p level, so 
that the spontaneous decay is neghgible. 

In the present scheme, ^3 and fJg correspond to iden- 
tical but counterpropagating lasers, and are thus inter- 
changeable. In the driving process from |1) to |2), an 
atom is equally likely to absorb a photon from either of 
the two lasers. Then, during the coupling from |3) to 
|4), it will preferably absorb a photon from the other 
laser, following by collectively decay. This is because for 
the competing process of absorbing a same-momentum 
photon, the phase matching condition for collectivity is 
not satisfied. Consequently, atoms driven by this chan- 
nel only decay spontaneously, the rate of which is a fac- 
tor 1/D <^ 1 smaller than that of the collective, phase- 
matching channel. 

The system's dynamics is solved by extending the rate 
equations ([l}-([6| to include intermediate pumping levels, 
such as levels of 6s6p,mp — —1/2,1/2 and 6s5d,mp — 
—3/2,3/2, as well as important side transition channels, 
including atom loss from |1) due to excitation by laser 
fli and CIq. The resulting rate equations are obtained 
in a straightforward manner very similar to the Cesium 
example derived explicitly in the appendix [B| 

We consider a cold (< 1 kelvin) spherical cloud of 
10^ atoms with diameter L = 26 /im, corresponding to 
an optical depth of Z? = 20. The six pumping lasers 
are chosen such that the Rabi frequencies and detunings 
of transitions denoted in Fig. |4] yield values of (all in 
units of GHz) = 5, ^2 ^ 16, = 2, = 10, 
= 35, fie = 2, Ai = 1000, A2 = 1200. With these pa- 
rameter choices, the effective Rabi strength resonantly 
coupling |1) to the level of \F — 3/2, nip — 3/2) is 
fioff = = 0.04GHz, obtained by adiabatically elim- 
inating the intermediate 6s6p level using fli,fl2 <C Ai. 
The same lasers couple |1) to another hyperfine level of 
|c) = \F = 5/2, mi? — 3/2), but blue detuned from res- 
onance by Afif — 4GHz. The effective Rabi frequency 
is riofi/2 <C ^hf- This, aided by that the Rabi fre- 
quency between |c) and |2) is 20,3 ~^ ^hf, suppresses 



the pumping from |1) to |c). This suppression necessar- 
ily provides a non- vanishing coupling channel from |1) 
to 1 2). Similarly, the coupling between |3) and |/) = 
\F = 5/2, —3/2), yielding an effective Rabi frequency of 
0.14GHz<C A/j/, is suppressed, too. We have numerically 
verified this analysis, where the population ratio between 
level \h) = \F = 3/2,mp = 3/2) and |c) is found to be 
4, while the ratio between |e) = \F = 3/2, mp = —3/2) 
and I/) is 14. 



(a) 




t fnsl ^ [ns] 



FIG. 4: (Color online) A realistic butterfly scheme using 
'^^^Yb atoms. Similarly to flgure|4]for Cesium atoms, figure 
(a) shows the level diagram, (b) plots the time evolutions of 
Npair and Nioss , and (c) draws the time correlation function. 



With the above observations, we have solved the ex- 
tended rate equations numerically, with results shown in 
figure [4] (b) and (c). In figure (b), the production rate 
of photon pairs is about 10^^ per second, while the atom 
loss from initial momentum modes is about lO^^s"^. The 
ratio of generated pairs Npair to lost atoms Nioss is then 
10. In figure (c), the time correlation function exhibits 
sharp peaks, showing strong pair correlation between sig- 
nal and idler photons. The bandwidth of signal and idler 
photons, measured by the reciprocal of correlation time, 
is about 15MHz, much smaller than the natural linewidth 
r4 — 0.2GHz. This is because the coupling between |3) 
and 1 4) is much weaker than the collective decay, so that 
the system dynamics is in the overdamped regime. 

The major atom loss is attributed to |1) atoms be- 
ing off-resonantly pumped to excited levels by fii and 
fli, and then decaying non-coUectively by emitting rogue 
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photons. These rogue photons are different in frequency 
from signal and idler photons by ~ 1000 GHz, which 
is much larger than the bandwidth of single and idler 
photons 0.02GHz). They can thus be easily filtered 
out from the photon-pair beams by optical frequency- 
selectors, and therefore will not contribute to the im- 
purity of collected biphoton beams. By excluding such 
rogue photons, the ratio of "good" paired photons to the 
remaining rogue photons, which are nearly frequency- 
degenerate with the photon pairs, turns out to be around 
20. For a twin beam of N pairs, this results in a fluc- 
tuation of in the number difference between the 
two beams. In contrast, the number-difference fluctu- 
ation for two uncorrelated coherent beams of the same 
size is y/2N . Hence, the generated twin beams in the 
present setup yield a considerable number-squeezing fac- 
tor of V2N/^/N/^ = 2^10. 



B. Cesium Atoms 

We now consider another implementation conflgured 
on the 852-nm line of D2 transition in Cesium atoms. 
The |2), |3) states are corresponding to \F = 3,nip = 
3), |F = 5,mF = 5), |F = 4,m_F = 4) levels, respectively. 
State |4) is consisted of degenerate \F = 4, mp = 4) and 
\F = 5, mp — 4) hyperflne levels. Driving from |1) to |2) 
is accomplished via a FWM process, through detuned 
intermediate levels of 6P3/2 and 605/2- The three driv- 
ing lasers, with Rabi frequencies fii, and fia, are tt, 
cr+ and cr_ polarized, respectively. In order to match 
the frequencies of signal and idler photons, the fis laser 
is blue-detuned by the ground hyperflne splitting of 9.2 
GHz. We note since the lifetime of the 6D5/2 level 1 
fis) is about 30 times longer than that of 6^3/2 30 ns), 
fl2 can be tuned near- or on- resonant without resulting 
in a faster atom loss than the photon-pair gain. The cou- 
pling between |3) and |4) is provided by a single resonant 
TT-laser. Signal photons are emitted as |2) atoms sponta- 
neously decay to |3), while idler photons are generated 
as |4) atoms collectively decay back to |1). Thus differ- 
ent from the toy level scheme plotted in Fig. [T] and the 
Ytterbium scheme shown in Fig. |4] the emitted photon 
pairs now have the same polarizations. Collecting these 
pairs directly generates a highly spin-squeezed beam in 
twin-Fock state, which is potentially useful for precision 
interferometry measurement at Heisenberg-limited sensi- 
tivity [36l [37] , as well as on-demand quantum teleporta- 
tion among single-atom qubits |38j . 

The system's dynamics is solved by extending the 
rate equations ([l])-([6]) to include all intermediate pump- 
ing levels, such as levels \a) = |6P3/2,mi? = 3) and 
\b) = \6D^/2,mp = 4), as shown in flgure|4|(a). Each 
level is associated with a corresponding spontaneous de- 
cay rate to account for background Rayleigh and/or Ra- 
man scattering of pumping lasers. Furthermore, there 
exist undesired side transitions due to cross-driving by 
pumping lasers. For example, atoms in |1) state are also 



driven by the laser to excited states, and |3) atoms are 
additionally coupled to |4) by the fii laser. These tran- 
sition are, however, far detuned from resonance. The net 
effects are then atom losses from relevant levels, and can 
thus be conveniently included in the rate-equation model 
by adding appropriate loss rates to corresponding levels. 
In this way, the multi-photon nature of driving and/or 
coupling is exactly treated via a dynamical model, while 
the background scatterings involving side transitions are 
calculated in a semi-exact manner, whose validity is jus- 
tifled by large detunings. We present the extended set of 
rate equations for the present Cesium scheme in appendix 

m 

To examine the performance of the Cesium scheme, we 
consider a spherical cloud of 10^ atoms with a diameter of 
L ~ AA /im, corresponding to D — 30. The temperature 
of the cloud is assumed to be well below 1 kelvin, so that 
Doppler broadening of level |2) is negligible compared to 
its natural line width. We choose the pumping parame- 
ters as (all in units of GHz): f^i ^2,^2 = 100, = 20, 
= 0.25, Ai = 300, A2 = 5, A3 = 9.2. We numerically 
solve the extended rate equations, with results shown in 
flgure|4](b) and (c). In flgure (b), the production rate of 
photon pairs is 1.7 x 10^"'^ per second, while the atom loss 
rate is 0.2 x lO^^s"^. The photon-pair gain to atom loss 
ratio is then ^ 9. In flgure (c), the time correlation func- 
tion exhibits sharp peaks, indicating strong pair correla- 
tion between signal and idler photons. The bandwidth of 
signal and idler photons, measured by the reciprocal of 
correlation time, is about lOOMHz, which is larger than 
the corresponding natural linewidth of 33 MHz. 

For the present setup, there are three dominating atom 
loss channels. The flrst is via spontaneous relaxation of 
atoms in the intermediate pumping level 1 61^5/2, "t-f = 
4), imparting emitting rogue photons at 917nm wave- 
length. The second is via background Rayleigh and/or 
Raman scattering of the fii laser by the populated |1) 
level, emitting photons at a frequency 300 GHz smaller 
than that of paired photons. The third channel is via 
spontaneous emission of |4) atoms to |3) and the ground 
\F — A,mp — 3) level, emitting rogue photons whose 
frequency is 9.2GHz smaller than the paired photons. 
All these frequency-differences are much larger than the 
bandwidth of generated photon pairs (~ O.lGHz), and 
can thus be flltered out by optical frequency-selectors. By 
excluding the aforementioned rogue photons, the ratio of 
"good" paired photons to the remaining rogue photons 
in the collected biphoton beams turns out to 25, corre- 
sponding to a considerable squeezing factor of 5\/2- 



IV. REABSORPTION ANALYSIS 

In any scheme using atomic ensembles for biphoton 
sources, photon pairs are generated inside an optically- 
thick atomic cloud. This raises the question of whether 
or not they are able to propagate out of the cloud with- 
out reabsorption. In some schemes |T21 |T51 [22 j sup- 
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FIG. 5: (Color online) A realistic butterfly scheme using Ce- 
sium atoms. Figure (a) shows the level diagram, (b) plots 
Npair and Nioss, and (c) shows the time correlation function 
for photon pairs emitted along ±S directions 



pression of this reabsorption occurs naturally, due to 
the existence of an inherent electromagnetically-induced 
transparency (EIT) window [51[51[3n]- For example, in 
Ref. |22i, where the generated photon pairs yield narrow, 
sub-natural bandwidths due to overdamped coupling, the 
idler photons are shown to propagate within the systems 
EIT window, thus suppressing the reabsorption. The sig- 
nal photons, on the other hand, are off-resonance with 
respect to excitation of the only populated |1) level, and 
thus in fact "see" an optically thin medium. 

In the Ytterbium scheme, similar EIT windows exist 
for both signal and idler photons as well. For signal pho- 
tons, the EIT window is formed by the |1) \a) <-> 



|^D2,r7i_F = 1/2) A-structure. The signal photons weakly 
drive the |1) \a) transition, while the resonant laser 
strongly drives the \a) ^ |^D2,m^ = 1/2) transition. 
Since i) the level (with a lifetime of 6700ns) is practi- 
cally metastable compared to ^Pi (with a lifetime of 5ns), 
and ii) the couplings of to other levels are negligibly 
weak compared to fis, the present A-structure is effec- 
tively mapped onto the three-level EIT model 15, \&j .22] . 
For idler photons, similar arguments apply, with the EIT 
window formed by the |1) |4) <-^- |e) A-structure. 
These, together with the fact that the generated pho- 
ton pairs yield a sub-natural bandwidth, give rise to sim- 
ilar EIT reabsorption-suppression effects demonstrated 
in previous experiments [iHl [22] , so that in the present 
setup the majority of signal and idler photons will be able 
to propagate out the atomic vapor. 

In the Cesium scheme, the signal and idler photons, 
which are both frequency and polarization degenerate, 
identically 'see' a EIT window formed by the |1) |4) ^ 
1 3) A-structure, where the transition |4) ^ |3) is reso- 
nantly driven by the strong laser, as in figure |5] (a). 
However, unlike the Ytterbium and previous schemes, 
here the bandwidth of photon pairs is of the order of 
1/lOns — O.lGHz, as from figure [5] (c). This bandwidth 
is broadened from the natural linewidth 0.03GHz. Thus, 
it is unlikely that the EIT-suppression effect still applies, 
as the bandwidth of photon pairs now matches the trans- 
parency window. On the other hand, the pulse duration 
of the photon pairs is {DTj^)~^, whereas any scattering 
event occurs at a rate not exceeding The proba- 
bility for photon scattering is then upper-bounded by 
(£)r4)-ir4 = l/Z? < 1. To verify this, in the follow- 
ing we use a simplified model to simulate the absorption 
of a single-photon pulse passing through a A-structure 
cloud. To account for the super-natural bandwidth of 
_Dr4, we approximate the single- photon Rabi frequency 
^y{t) with a "toy" envelope function of 



ny{t) = f7oe"-°^''*/^sin(r24</2) 



where the coefficient 



8IVD 



Ntt 



(24) 



(25) 



is determined such that the underlying EM field yields 
the single-photon energy /ic|K|. The atomic dynamics of 
this A system is governed by a set of rate equations. 



dt 
_d 

dt 

d 

~dt 



-Na 



(ftygii - c.c) + T^iPNi{Ni + 1), 
{^4943 - c.c) -I- TispN^Na + 1), 

^(^^4^^43 + i^yQAl - C.c) - (r4 + ^(r4liVi + r43iV3))^4, 



(26) 
(27) 
(28) 
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dt 
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dt 



en = 



Sis 



631 



N3) 

.^4 



' 2 



43 



Ni) 



N3) 



r 



41 



^A^4 ^31, 



£43 

2 

£41 
2 



PN3 



^43, 



(29) 
(30) 
(31) 



obtained in a similar manner with Eq. ([T|)-(|6|. Here, Ni, 
N3 and A^4 are populations at level |1), |3) and |4). gij 

= 1,3,4) is the coherence between level \i) and 
Tij is the spontaneous emission rate from \i) to \j) and 
Ti = Tij is the natural linewidth of level \i). Starting 
with all atoms in the |1) state, the above equations of 
motion are solved numerically. Because photon losses 
correspond to scattering atoms out of the A system, the 
loss probability Pioss of the photon passing through the 
medium is given by the atomic population reduction. 



,(T) =iVi 



t=o 



{Ni + N2 + N3) 



t=T 



(32) 



A plot of Pioss{T) is shown in figure [6] for parameters 
used in section [riIB| We find the ultimate loss probability 
Pioss{oo) K, 0.02 for both signal and idler photons. This 
result is consistent with the upper- limit of 1/_D w 0.03 as 
from the time-scale argument. In practice, this photon 
loss due to absorption will lead to a 2% degradation of 
the pair correlation. 




FIG. 6: The loss probability Pioss{T) as a function of time T 
for both signal and idler photons in the cesium scheme. 



V. CONCLUSION AND DISCUSSIONS 

In this paper we have proposed an omnidirectional 
biphoton source based on collective emission in an atomic 
vapor. This is accomplished by employing multi-photon 
excitation and 'Doppler-free' pumping. Our scheme ben- 
efits from the elimination of background Rayleigh scat- 



tering via dipole-selection rules, and the fact that phase- 
matching can be fuUfilled in the full 47r solid angle. We 
have demonstrated an achievable photon-pair brightness 
of several orders of magnitude greater than the best re- 
ported results. Our scheme has the unique feature of 
strongly suppressed atomic re-thermalization, thus al- 
lowing implementations in 'hot' vapors as well as ul- 
tracold samples such as BEC's, performed nondestruc- 
tively. This may lead to nondestructive in situ imaging 
of condensates, with an anticipated atom-counting preci- 
sion below the standard quantum limit. This is because 
potentially, more photons than atoms can be generated 
before destroying the condensate. Furthermore, we found 
strong time and polarization correlation between signal 
and idler photons. We also provided two realistic im- 
plementation of the present scheme, using Cesium and 
Ytterbium atoms. In both schemes, the generated pho- 
ton pairs can propagate through the optically thick cloud 
without being reabsorbed. 

In the present scheme, the impurity of collected bipho- 
ton beams are attributed to i) non-collective emissions 
from level |4), and ii) background Raman /Rayleigh scat- 
tering of pumping lasers by the |1) atoms. The branch 
ratio of collective to non-collective emission is D, i.e., 
the optical depth of the atomic sample. One then might 
seek to improve the beam purity by employing samples of 
higher optical depth. The difficulty is that this requires a 
stronger coupling between |3) and |4), in order to main- 
tain a strong time correlation between signal and idler 
photons. In both Yb and Cs schemes, this will unavoid- 
ably lead to a stronger background scattering, adding to 
the beam impurity. 

Hence, an optimal optical depth is in fact obtained 
by balancing the two impurity sources of non-collective 
emission from level |4) and Raman/Ray Icigh scattering 
of pumping lasers by the |1) atoms. To strongly drive 
the 1 3) ^ 1 4) transition while avoiding significant back- 
ground scattering, it is clearly necessary to seek appro- 
priate level structures where the coupling is resonant for 
the |3) ^ |4) transition, but far detuned from exciting 
|1) atoms. An example is seen in the Cesium scheme 
in section [ill B where the coupling laser Q.^, while reso- 
nant with the 3) <-*■ |4) transition, is red-detuned from 
the |1) excitation due to ground hyperfine splitting. For 
further improvements, one may seek level scheme con- 
figured on the fine structure to exploit the extremely 
large level splitting (> THz). Feasible atomic species in- 
clude those having no ground state hyperfine splittings. 
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such as ^^°Sn and ^^Si. For these atoms, we have found 
butterfly level structures similar to the cesium scheme 
shown in figure [S] (a) , but configured on fine structure 
levels, where the spacing between level |1) and |3) is now 
10"* ~ lO'^ GHz, instead of 9.2 GHz as for Cesium. These 
schemes are expected to have very low background scat- 
tering. The disadvantage, however, is that because the 
omnidirectional phase-matching condition relies on pro- 
ducing (nearly) identical-wavelength signal and idler pho- 
tons, the 1 2) level must be blue-detuned by an amount 
of ground fine splitting of lO''^ ^ 10^ GHz, which would 
then require ultra-intense driving lasers, making it diffi- 
cult implement a cw scheme experimentally. Future work 
will therefore study the feasibility and applicability of a 
pulsed driving scheme. 

If desired, level structures similar to the Cesium 
scheme in figure |5] (a) can be applied to other Alkaline 
atoms, such as Rb and Na. However, they have smaller 
ground hyperfine splittings than Cs, and will thus be sub- 
ject to stronger background scattering. To overcome this 
difficulty, one possibility would be to increase the spht- 
ting by applying a static magnetic field, making use of the 
fact that level |1) and |3) are exposed to different Zeeman 
shifts. For a rough estimation, increasing the splitting to 
10 GHz only requires a magnetic field of ^ 0.1 T, sug- 
gesting experimental feasibility. 

Lastly, we note that as the atom loss fraction becomes 
significant, the system will enter a more complicated 
regime, where the system tries to equilibrate, potentially 
resulting in macroscopic occupation of all ground hyper- 
fine sub-levels, and presumably strongly diminished pair 
correlations. If this is the case, only a small number 
of photon pairs will be generated, but in a very short 
bright initial burst. Thus, the butterfly system may be an 
excellent source for generating highly number-difference 
squeezed twin pulses for quantum interferometry. Fur- 
ther studies of the butterffy system in this equilibrium 
regime remain as a future task. 

Acknowledgement: we thank A. Leanhardt for helpful 
discussions. This work is supported in part by National 
Science Foundation Grant No. PHY0653373. 



APPENDIX A: DERIVATION OF RATE 
EQUATIONS 

In this appendix, we describe how one can derive the 
rate equations ([l])-(|6]) from first principles, based on a 
Markovian treatment of the optical field. We begin by 
deriving a coarse-grained version of the Heisenberg equa- 
tion of motion for an arbitrary operator, A. This ap- 



proach is equivalent to the standard master equation ap- 
proach, but eliminates the density operator, as it is an 
unnecessary intermediate step when the goal is to derive 
equations for expectation values of observables. 

The standard quantum time propagator is given by 
U{t) — e^*^*. Defining in the usual way, the Heisenberg 
picture operator, A{t,tQ) ~ W{t — to)AU{t — to), it the 
follows that A{t + T, to) ^U\t + T- to)AU{t + T-to), 
so that 



dt 



A{t,to) 



lim - ( JjUt- 

T^O T V 



to)AU{t- 



U\t~to)AU{t~to)) 



We can then set to — t to arrive at 



W{t)AU{t)- A 
dt ^ T 



(Al) 



(A2) 



where A is now the Schrodinger picture operator, and t 
must be chosen sufficiently small. Taking the expectation 
value of this equation with respect to the state of the 
system at time t, will then yield the time-derivative of the 
expectation value at time t. The Markov approximation 
then consists of letting t be small compared to the system 
evolution timescale, but long compared to the dephasing 
time of the reservoir. 

Using this approach, we now derive a generic equation 
of motion for the toy level scheme of Sec. |II A| We start 
from a generic system-reservoir model. 



(A3) 



where the reservoir consists of a bath of bosonic field 
modes, governed by a Hamiltonian of the form 



(A4) 



where annihilates a bath particle with wave- vector, k, 
and frequency, Wk, and LOg is the resonance frequency for 
the system-reservoir interaction. We assume an interac- 
tion operator of the form 



E 



(A5) 



where Ck is an unspecifi ed s vstem operator. 

Now let us evaluate (A2| for the case A ^ S, where 
5 is a system operator only. To second-order in Vsr, this 
gives 



dt T 



WoiT)SUoiT) - S 



dto 



dtiC/T(il)V;.(ti)t/J(T - ti)SUoiT - t2)K,(t2)l/o(<2) 
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1 

T JO 
1 
T 



dt2 / dhUoir - t2)SUo{T - <2)T4,(t2)C/o(t2 - h)Vsr{h)Uo{ti) 

Jo 

dt2 f ' dhU^oih)VsAh)U^oit2 - tl)trih)Ul,iT - t2)SUo{T) 

Jo 

I 



(A6) 



where Uo{t) ~ g-i{Hs+Hr}t ^ g^^^ have neglected the 
first-order terms as they will vanish when we trace over 
the reservoir degrees of freedom. For the Markov approx- 
imation, we first take 



C/o(t) « (1 - iTH,)e- 



(A7) 



and then perform the reservoir trace. The time integrals 
are then handled via, 



(A8) 



which gives 
d 



dt 



S 



V 

k 

i 



Ck + 4 



S,cu 



[6.. 



s 



(A9) 



If one desires, the equation of motion for the system den- 
sity operator, ps, can be obtained as a special case of 
(|A9|) , with the substitution ^ 



4^ , as the state of the 

at ' 



system must be evolved backwards in time to obtain the 
state in the Heisenberg picture. 

For the toy model depicted in Fig. [l] the system 
Hamiltonian is 



N r 

E 

2 



iK-f j 



|2)(l|0j 



-i'K.-rj 



|l)(2|0j) 



|4)(3|0j. 



|3)(4|0j) 



, (AlO) 



where f ^ is the position operator of the j*'* atom, and 
the state \m)0j indicates that the j*'* atom is in inter- 
nal state |m). The system Hamiltonian (AlOl is given 



for a frame rotating at the system resonance frequency, 
u!s, corresponding to the frequency of the |4) |1) and 
1 2) 1 3) transitions. The system-reservoir interaction is 
described by the system operators 



Ck 



N 



(|l)(4|0j + |3)(2|0j) 



(All) 



For the initial state of the system, we assume that each 
atom is an internal state |1), and occupies a single box 
eigenstate, q, so that 



N 

n i^j)^-?' 



|i)0i- 



(A12) 



The initial momentum of the j*'' atom, is to be cho- 
sen at random from the Boltzman distribution. For large 
enough N, this will reproduce the results of thermal av- 
eraging in a single realization. We then introduce the set 
of states 



q)0j <S) \l)0j 
K + q)0j(g>\2)0j, 
K + q)0j0 |3)0j, 
q)0j ® |4)0i, 



(A13) 
(A14) 
(A15) 
(A16) 



\Uqi)0J iQj 

\Uq2)0j hj 
\Uq,3)0j ■■= h] 

\uq4)0j := |qj 

so that for the j*'* atom, |moi)0j) and |wo2)0j are the 
initial state and the state after interacting with the drive 
lasers, respectively. The states |itq3)0j and |uq4)0j, then 
correspond to the state after emission of a signal photon 
with momentum k w — q, and the state after then ab- 
sorbing a sequence of photons from the coupling lasers, 
respectively. Emission of an idler photon with k w q, will 
then return the atom to the |uoi)0j; whereas emission of 
a rogue photon will transform the state into |uqi) with 
q^O. 



With respect to the states (Al3l-(A16), the system 
Hamiltonian becomes 



E 



2 



qlq2 



s. 



t 

qlq2 



s, 



q3q4 



s. 



t 

q3q4 



(A17) 

where we have introduced the generic atomic transition 
operators. 



Sn 



N 

E 



u^,){uy\0j, 



(A18) 



where ^, e {{ql}, {q2}, {q3}, {q4}} are composite in- 
dices. Similarly, the interaction operators take the form 



Ck 



where 



E/(k- 



/(k+q-Q) 



Q) ^, 



'qlQ4 



q3Q2 



(A19) 



1 

V 



-ikr , 



IQ> 



^3^g-»(k+q-Q).r (^20) 



is the static structure function of the sample. 

The operators we are interested in are all one-body 
operators of the form (Al8l. In order to evaluate (A9l 



for these operators, we will need to evaluate commutators 
of the form. 



M/3 



(A21) 
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Taking the expectation value of the equation of motion 
(A9) will then require us to compute the expectation 
value of the bilinear operator (S^^Safs)- Our strategy 
for dealing with these terms will make use of the under- 
lying exchange symmetry which leads to the emergence 
of collectivity in the emission properties of the sample. 
We note that the initial state (A12l is not symmetric 



under particle label exchange, due to the dependence of 
the state |woi)0j on q^. However, as the set of states 
{|"qi)0i, I"q2)0j, {\uq3)0j}, {|uq4)0j}} forms a complctc 
basis, all expectation values of operators of the the form 
(AI8I will depend only on inner-products of the form 
{u^\u„)0j, which are in fact independent of the q/s. For 
example, with m,M & {1, 2, 3,4}, we have 

(wqm|tiQM)0j = ("T-|0j ® (qi+q|0j 1^)0^ ® |qj+Q)0j 
= W/(q-Q) (A22) 

which does not depend on q^. This means that we can 
make any choice we like for the set of q^ 's, without affect- 
ing the rate equation dynamics. We note that this occurs 
in part because we have neglected the kinetic energy of 
the atomic center-of-mass motion in our system Hamilto- 
nian ( AlO I, an approximation valid when the lifetime of a 



collective excitation is short compared to 1 / {u d{T)) , 
where lur = fiK^/2M is the recoil frequency, and Yd is 
the Doppler line- width ( 22 1 . 



One option is then to set all the q^-'s to zero, in which 
case our initial state becomes explicitly symmetric under 
exchange of particle label exchange. Equivalently, one 
can recognize that the q^'s are redundant with the par- 
ticle labels, so that exchange of the q^s is part of the 
underlying symmetry. As the Hamiltonian is also sym- 



metric under this form of particle exchange, it follows 
that the state of the full system will remain symmetric 
as it evolves in time. The full state of the system + 
reservoir can then be written as 

\ipSR{t)) = ^ci^ii,. . . ,fiN;r)\Uf,^)0l<»\Uf,^)02 



/ii,...,/jjv r 

®|u^3)03 . . . (g) |Mp„)0iV ® \r)0R, (A23) 

where r is a composite index which sums over all states 
of the reservoir, and the state \r)0R lives in the reservoir 
Hilbert space. Particle exchange symmetry then requires 
that c(/ii, . . . , /ijv; r) be invariant under exchange of any 
two /ij's. 



With the symmetric state ( A23 1 , we can now evaluate 



the expectation value of the product of two operators of 



the form (AI8I 



N 



{S^^Sc^fi) ^ d^c^iSf^p) + Y {S^,^0jSo,p0j), (A24) 



S3, j=i 



where 

Si,u0j \u^){u^\0j, (A25) 
and with ^ — qm and v — QM, we have introduced 



Oq,QOm,M- 



(A26) 



To illustrate an important consequence of exchange sym- 
metry on the bilinear terms in ( A24 1 , we consider first 
the j=l, J=2 case, 



(^^,0l5„^02) = (^/.(t)IV0l^a/302|V'(i)) 



/i3,...,/xjv r 

= {S^p0lS^,02), (A27) 

I 



5^i,v to elimi- 



where we have used the relation {Ufj_\ui, 
nate these inner products. Because the state is symmet- 
ric under exchange of particle labels, Eq. ( A27[ ) can be 
generalized to {Sf^y0jSap0J) = {S^f30jSau0J), for j =^ J, 



which gives us 

{Sf,Ja0)=d^,a{Sf,p)+ Y {S^,fi0jSa^0j), (A28) 



N 



Bj, J = l 



as an equivalent alternative to (A24l. This result will 



allow us to express the equations of motion for number 
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operators in terms of products of number operators, as 
opposed to products of coherence operators. 

The equation of motion for the expectation value of a 
system operator of the form (AI8I is then 



A. 
dt 



(4 [ 



Ck) 



(A29) 



We can evaluate the commutators via Eqs. ( |AT7| , ( |Al9l ) 
and (A2I1, resulting in 



k' M^" 



(A30) 



^/(k+q-Q) (5 

(A31) 



q,Q 



and 



2 



' 2 



q 



q2'J/^ql 



q 

+ <5/i,q3'S'q4l/ — ^l/,q4*S';^q3 



(A32) 



At this point, the variables used in the rate equations 
pT|6|, can be precisely defined as 



Ni 




( om ni ) . 
\ u lu 1 / 1 


(A33) 


N2 


= (^02)- 


('S'0202), 


(A34) 


Q12 


— {80102), 




(A35) 


iVq3 


= (^q3> = 


('S'q3q3), 


(A36) 


iVq4 


= (^q4) = 


('S'q4q4), 


(A37) 


£'q34 


= (5'q3q4)- 




(A38) 



The occupation numbers A'^i and A'2 count the number 
of atoms in internal states |1) and |2), respectively, that 
have not been 'lost' via emission of a rogue idler pho- 
ton. The occupation numbers iVqs and A^q4 count the 
number of atoms in internal states |3) and |4), which are 
displaced in momentum space by K + q and q, respec- 
tively, relative to their initial momenta, corresponding 
to their having emitted a signal photon with momentum 
k « — q. The coherence terms Q12 and gq34 describe the 
coherence generated by the driving and coupling lasers, 
respectively. 

Beginning with A^i, we can derive its equation of mo- 
tion from ( A29 1 by setting ^ = u = 01. For the commu- 
tators (|A30rfA32l we find 



cj^, S'oioi 


= ^r(k-q)5q40l, 
q 


'S'oioi, Ck 


= ^f{^-q)Soi,ii, 




q 


Hs, S'oioi 


— (^'5oi02 ~ 'S'0201^ 



(A39) 
(A40) 
(A41) 



which leads to 



-TVi ^Y.^8(uj^-UJs)Wf /*(k-q)/(k + Q'-Q) ((5q401^Q'lQ4) + (^q40l5Q.3Q2; 

k q,Q,Q' 

I 



i ^ £'12 + C.C. (A42) 



With the help of (A27l we see that 



N 

('S'q401'5'Q'lQ4) = '^Q',o('5'q4Q4)+ ^ ('S'q4Q40j<5'Q'lO10>^) • 

flj, j=i 

(A43) 

To implement the approximation that all atoms that emit 
rogue idler photons are permanently 'lost', we simply 
evaluate expectation values under the assumption that 
there are no atoms in the states |uqi) and |uq2) for q 7^ 0. 
This allows us to make the simplifications 



('S'q40l5'Q'lQ4) — <^Q\0 (Afq4 + (iVq4iVoi; 



(A44) 



and 

('S'q40l5'Q'3Q2) = '5Q,o(5'q401<5'Q'302)- (A45) 

Inserting these into ( |A43 ) and making the approximation 
r (k - q)/(k - Q) « |/(k - q)|2 ^Q,q , (A46) 
then gives 
d 



dt 



J2 ((^q4^01> + A^q4 + (S'q40l5'-q302>) 



i-z-gi2 + C.C. , 



(A47) 
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where we have introduced the branching ratio 



/3q 



^7r<5(a;k-c..)|5kn/(k-q)P, (A48) 



Turning now to the approximation (b), we begin by 
introducing the reduced two-body density operator 



(A57) 



where 



27r(5(wk 



k 



I5k| 



(A49) 



is the spontaneous emission rate for states |2) an d |4) , 
which we have impUcitly set equal by making Ck (All) 
symmetric with respect to the two transitions. 

The next step is to show that (a) (NqiNoi) w 

(/Vq4)(7Voi), and (b) that (S'q4oi5'_q302) < (iVq4^oi), so 
that it can be safely neglected. In order to verify (a) and 
(b), it is useful to define the excitation- number operators 



and 



Nd := iVq3 -I- Nq4- 



iVe:=5]iVq 



(A50) 



(A51) 



with the number of non-excited atoms then given by 

No:= Ni+N2^ N - (A52) 

In terms of probabilities, the unproven approximation (a) 
can be re-expressed as 



NiN - l)P2(01,q4) « iV2Pi(01)Pi(q4), 



(A53) 



where P2{n,i') is the joint probability that for any or- 
dered pair of atoms, the first one will be in state 
and the second in state \ui,), whereas Pi(/i) = N^^/N is 
the bare probability that any given atom will be in state 
\Ufj,). According to Bayes theorem, we have 



P2(q4,01) =Pi(01|q4)Pi(q4) 



(A54) 



where Pi(01|q4) is the conditional probability to find a 
particular atom in state |moi) given that another particu- 
lar atom is in state |uq4). Knowing that the second atom 
is in state |Mq4) means that of the remaining iV — 1 atoms, 
the average number of excited atoms is now Ne ~ I, or 
equivalently, out of the remaining iV — 1 atoms, the av- 
erage number in state |uoi) is still TVi, Thus we see that 



Pi(01|q4) 



N- 1 
This leads to the result 



Pi(Ol) 



N 



N - 1 



(A55) 



(A56) 



i.e. the factorization (a) is exact. We note that this would 
not be obtained without assuming a symmetrized wave- 
function with fixed total atom number, which allowed us 



to replace (a) with Eq. (A53l 



so that (S'q4oi<5'_q302) = P2(q4, -q3;01,02). If we as- 
sume that the probability to have iVq > 1 is negligible, 
it follows that the matrix element p2(q4, — q3; 01, 02) is 
a measure of the coherence between the A'q = iV_q = 
manifold, and the A'q — iV_q = 1 manifold. 

To understand the origin of such coherence in the sys- 
tem, let us start from the initial state = Il^i Koi), 
and assume that at time t — 0, the driving and coupling 
laser beams are turned on. The atoms will adiabatically 
follow the ground state of the dressed system, so that 
immediately after the fields are turned on, the state of 
the system will be 



N 



Ho) = n 



(A58) 



This dressed state will decay by emitting a signal photons 
at the rate T^ff — Nil'^/T. The signal photons are dis- 
tributed over the many recoil modes of the system, with 
branching ratio (3 per mode. Thus from the perspective 
of a single iVq manifold, the time delay between signal 
photons is Tq = (3 /Terr = r/(f^d-D)i where we have used 
D « Np. 

Let us assume that the first photon is emitted along 
k « — q, causing the system to jump from the iVq — 
A^_q = manifold to the iVq — 1, A^-q = manifold. 
The normalized state immediately after this quantum 
jump is given by 



C-qlV'o) 



(Vo|cLqC_q|V'o)l/^ 



(A59) 



By neglecting rogue-photon emission and making the ap- 
proximation (A46l, we can simplify (A19l to 



Cq — 5'oiq4 + S- 



q302, 



which leads to 



IV'l) = -«^^'5q3O2|V'0)- 



(A60) 



(A61) 



This state will live for time ti ^ — , after which a 
second photon will be emitted along the k « q direction. 
With t — corresponding to the emission of the second 
photon, the state of the system at time t later is then 



\Mt)) = 



where the propagator 



Uc{t) = exp 



nrt 



q3q4 



q4q3 



(A62) 



(A63) 
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describes the Rabi oscillations between states |3) and |4), 
that occur when ilc > TD/2. 

We can take Uc{ti)Sq302 ~ Sq402, as Uc{t) mixes levels 
1 3) and |4), and the highest probability of emission occurs 
when the excited atom is in level |4). This leads to the 
result 



r 



5*0102 + C^c(^)'S'_q302'S'q402 ) IV'l 



(A64) 

which shows that the state |V'2(^)) is in fact a coherent su- 



perposition of a state with iVq 



0, corresponding 



to the second photon being an idler photon (emitted on 
the |4) — > |1) transition) and a state with iVq = N^q = 1, 
corresponding to the second photon being a new signal 
photon (emitted on the |2) — > |3) transition). The life- 
time of this coherent superposition state is Tcoh ^ T^, 
the timescale on which either (i) a second pair of pho- 
tons will be emitted, confirming that the second photon 
was a signal photon, or (ii) no additional photons will be 
emitted, confirming that the second photon was the idler 
photon. In either case, the system will collapse back onto 
a state with A'q — iV_q = 0. 

From this analysis, we now see that 



('S'q40l5'_q302) 



To 



( V'2 (i) I <S'q401 5'-q302 I ^"2 (*) ) 



(A65) 



where the factor T^/Tq = ri^/F^ is the probability to find 
the system in state |V'2(i))- To co mpar e this to the term 
(7Vq47Voi) = Ar2Pi(01)Pi(q4) in ( |A47| , we need to esti- 
mate the single particle probabilities, -Pi(/i). Estimates 
for these probabilites can be found from the equilibrium 
condition Rg — Rj. For ^ N and ilc > FZ3, we can 
make the simple estimates Rg = N2T and Rj = ^NeTD. 

With N2 ^ ^Ni and = 1 ~ Nq, we find A^i = 



and A^qs « A'q4 



N2 



r2+a2(i+2/_D)' 
^PNe, which gives us 



Pi (01) « 1 - 



p2 ' 



{2/D)niN 

r^+n^ji+2/D)' 



(A66) 



and 



Pi(q4) 



Thus we see that 



(A^q4A^01 



^1 

F2 N' 



f]2 

— N 

p2 ' 



(A67) 



(A68) 



which is larger than (S'q4oi5'_q3oi) by a factor (F/fi^)^. 
For th e par ameters used in the numerical simulations of 
Sec. (II B I, this is a factor of 100. Keeping only the 
dominant term, (NqiNoi), then leads to 



dt 



Al -F^/3qA^q4(A^l-|-l) 



«^ {Q12 



C.C.), (A69) 



which for F4 = F and /3q4 = /3q, validates the rate equa- 
tion ([T|). Using the general methods outlined in this ap- 
pendix, the remaining rate equations ([2|-(|6| can be de- 
rived as well. Clearly deriving these rate equations from 
first principles is highly non-trivial. However, once their 
validity is established, the fact that they follow an es- 
tablished form allows one to write them down directly, 
rather than re-derive them for each particular model. 



APPENDIX B: EXTENDED RATE EQUATIONS 

In this section, we give the extended rate equations 
we have used to model the Cesium scheme presented 

to |2) 



m section 



IIIB 



In the scheme, driving from |1) 
is accomplished via a three-photon process, through de- 
tuned intermediate levels of \a) = |6P3/2,™_F = 3) and 
\b) = |6£'5/2,TO_F = 4), as shown in figure Is] (a). With the 
conventions of Nj being atom number in sTate | j) and Qjh 
being the coherence between \j) and \h), the equations 
are obtained as 



A^i 

Na 



dt 

d 

dt 
d 

~dt 



- ^ {^iQai - c.c) - FiA^i + a4ir:^^ /3k4A^k4(A^i + 1), 

k 

= {^iQal - C.c) + ^ {^2Qba " C.c) - F°A^a, 

= {^2Qba - C.c) - ^ {^iQb2 - C.c) - F°A^f„ 

= J (^^3^62- C.c) -F2^/3k2^2(^k3 + l), 



(Bl) 

(B2) 
(B3) 
(B4) 
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dt 

d 
di 

A. 
di 

di 

d 
dt 

d 
di 



Qal 



Qba 



Qbi 



Qb2 



Qa2 



621 



dt 



^zAi-^(r[J + ri) + ^a4ir°X^/3k4Afk4j gai 



n 



^^{Na - m) - z^Bbi + ^^Ql2 ~ *(A2 + Ai) + -(r» + r"') q^. 



^2 ^1 ^3 

.f^3, 



1+ ^^a4irJl^/3k4A^k4- ^(Fi+rg)^ f-bi+zAaPfci, 



«^(iV2 - iVb) + + |^i(A3 - A2) - - \^2 /3k2(iVk3 + 1) j gb2, 



2 



^^2 ^3 
i^Pl2+l^Qb2 - '>'^QI 



.^1 



-^^Pa2+^^Qbl + 



\a + (^*(Al + A3) - - X^r2/3k2(^k3 + 1)^ Qa2 

I A3 - - ^ (r2/3k2(A^k3 + 1) - a4ir:i/3k4iVk4 



2 (f^4ek43 - C.C) + r2/3k2^2(A^k3 + 1) - r3Ar3, 
(f^4£'k43 - C.c) - r°iVk4(a4l/9k4A^l + 1), 

«^(A^k3 - A^k4) + ^ (^r2Aik2A^2 - r3 - r!^(a4iAik4A^i + 1)^ ek43, 



(B5) 
(B6) 
(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 
(B12) 
(B13) 



assuming all fi's are real. Here, Aj {j = 1,2,3) is the 
detuning indicated in figure [s] (a). r°, r|j, r2 and 
are the natural linewidths of states \a), \b), |2) and |4), 
respectively, a^i is the branch percentage for |4) to spon- 
taneously decay to |1), so that r4 = a4ir4 is the sponta- 

neous decay rate from |4) to |1). Fi = ^j^r^ is the effec- 
tive loss rate of |1) atoms, due to being excited by the de- 



tuned laser fl^ to the unstable upper level |a). Similarly, 
r3 = so Ta is the effective loss rate for atoms in 



4(Ai-A3)2-"4 

1 3), induced by laser fii. We note the above rate-equation 
model has only incorporated dominant atom loss mech- 
anism for each atomic level. 
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